Strong Coupling of a Cavity QED Architecture for a Current-biased Flux Qubit 
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We propose a scheme for a cavity quantum electrodynamics (QED) architecture for a current- 
biased superconducting flux qubit with three Josephson junctions. The qubit operation is performed 
by using a bias current coming from the current mode of the circuit resonator. If the phase differences 
of junctions are to be coupled with the bias current, the Josephson junctions should be arranged 
in an asymmetric way in the qubit loop. Our QED scheme provides a strong coupling between the 
Hux qubit and the transmission line resonator of the circuit. 
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I. INTRODUCTION 

Superconducting qubits coupled to a quantum mechan- 
ical harmonic oscillator have reached the strong coupling 
region, where the coupling strength is much larger than 
the decay rates of the cavity and the qubit. The cir- 
cuit quantum electrodynamics (QED) scheme has been 
applied to the superconducting charge qubilii'^, phase 
qubit-3^, and flux qubit *"^. Recently, in Refs.^"^ the phase- 
biasing scheme was studied to provide a strong coupling 
strength between the resonator and the flux qubit by 
sharing the qubit loop with the resonator. This galvanic 
coupling, however, has difficulty in switching on and off 
the coupling between the qubit and the resonator, which 
is essential for a scalable design. In this study, by in- 
troducing a current-biasing scheme for the flux qubit, 
we offer a new circuit QED architecture, where the flux 
qubit is coupled with the current mode of the resonator, 
to obtain a strong coupling between the flux qubit and 
the resonator without the galvanic link. In this scheme, 
the flux qubit is biased by the oscillating current mode 
of the resonator. 

The current-biased dc-SQUID qubit (phase qubit )^"— 
is controlled by an electric field, which provides a fast 
operation. On the other hand, the flux qubili^iii^ is oper- 
ated at the optimal point where the first-order phase fluc- 
tuations vanish. The present current-biased flux qubit 
operation is also performed at a point optimally biased 
with respect to both the bias current and the external 
magnetic flux, which may provide a long coherence time. 
An oscillating bias current induces a Rabi oscillation be- 
tween qubit states. If the bias current operation of the 
flux qubit is to be performed, the number of the Joseph- 
son junctions in the flux qubit loop should be three (in 
general, an odd number), and the junctions should be 
arranged in an asymmetrical way. 



II. CURRENT-BIASED FLUX QUBIT 

Figure [Tfa) shows the current-biased flux qubit, where 
the current flows through the three- Josephson-junction 
qubit loop. In the circuit QED architecture with the flux 



qubit, which we will discuss later, the bias current comes 
from the oscillating current mode of the resonator. For 
the time being, we consider an externally applied oscillat- 
ing bias current. The Hamiltonian of this system can be 
derived semiclassically by using the quantum Kirchhoff 
relation. 

First of all, let's consider a superconducting loop with- 
out a Josephson junction, where the usual fluxoid quan- 
tization condition readsi^ 



- $t + {nic/qc) 



■ dl = 7i<I>o, 



(1) 



with Vc being the average velocity of Cooper pairs, qc 
being 2e, and TOc being 2TOe- The total magnetic flux 
threading the loop <I>t is the sum of the external and the 
induced flux <I>t = $cxt + *&ind with the superconducting 
unit flux quantum $0 = h/2e. Then, the condition of 
Eq. ^ is written as 



kl = 27r{n + ft), 



(2) 



where ft = $t/$o = / + find with / = $oxt/*o and 
find = $ind/*&Oi I IS the circumferencc of the loop, and k 
is the wave vector of the Cooper pair wavefunction. 

For the current-biased flux qubit in Fig. [ija), the flux- 
oid quantization condition in Eq. ([T]) is changed due 
to the phase differences 0i's in the circumference of the 
qubit loop as follows: 

{ki+k2)-=2Tr{n + ft)-^i-<j>2-(f>3- (3) 

Usually, for flux qubits, the contribution (fci -f k2)l/2 
is negligible in Eq. ([3]) which then can be reduced to 
2Tr{n-\-ft)—(j)i—(j)2—4'3 = 0- Rigorously, however, we keep 
this term for the time being to derive the Hamiltonian of 
our qubit system. 

The induced flux $ind is given by <i>ind — Ls{Ii — /2)/2 
with the self inductance L^, and the current 1^2) in Fig. 
(Ha) is 



h{2) = T{'ncAqc/mc)hki(2), 



(4) 



with He being the Cooper pair density and A the cross 
section of the superconducting loop. Then, by using 



Eq. dH) and *o = h/q^, find = *ind/$o can be rep- 
resented as /ind — —{Ls/LK){l/'2){ki + fc2)/27r with 
Lk = mJ/Aricqc being the kinetic inductanceJ^. With 
this relation the condition in Eq. ^ becomes 
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^^ -(5) 



The dynamics of the Josephson junction is described 
in the capacitively-shunted model, where the current re- 
lation is given by / = —IcSiiKp + CV with the capaci- 
tance of junction C, V — dV/dt, and voltage V across 
the junction. This relation can be rewritten by using the 
Josephson voltage-phase relation V — — ($o/27r)0 as 

- {ncAqc/mc)hki = -I„8m(j)i - Ci($o/27r)(?!)i, (6) 

with the critical current of Josephson junction Id — 
27r_Eji/$o {i = I7 2, 3) and the Josephson coupling energy 
Eji. Then, from Eqs. ([5]) and (jS]) with ko — ki — k2, we 
have the quantum Kirchhoff relation 
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Cichi = 
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47r 



(7) 



where /q — —{ncAqc/mc)hko. Here, for i — 1,3, the last 
term of Eq. ([7]) is -|-$o^o/47r whereas for i = 2 the sign is 
reversed as — $o/o/47r. For usual flux qubits, Lk/L^ ^ 
0.01—; thus, hereafter Lk is neglected for simplicity. 

The equation of motion, Eq. ([7]), can be obtained from 
the Lagrange equation {d/dt)dC/d(l)i — dC/d(l)i = with 
the Lagrangian 



^{4>i,4>i 



E- 



c. 



^1 - C/eff({0^}), 



(8) 



Ues{{4>i\) = ^Ej,{l-COS(j),) 
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(9) 



where the first term in Eq. ([5]) comes from the charging 
energy of the qubit system, Ec = J2i=i Ql/'^^i '^ith 
Qj — Ci($o/27i')0i. The second term of Eq. (0) has a 
finite value owing to the asymmetry of the qubit loop, 
giving rise to the coupling between the bias current and 
the flux qubit. 

In the usual experiments for fiux qubits, $q/2Ls ~ 
0{lQ^Ej) is much larger than the other energy scale 
in the Lagrangian of Eq. ([8]). Hence, the last term in 
Eq. ^ can be removed, leaving the constraint g{(j)i) — 
01 + 02 + 03 ~ 27r(n -|- /) = 0, which is the familiar fiux 
quantization condition. The Lagrangian then has the ef- 
fective potential 
3 
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FIG. 1: (Color online) (a) Current-biased flux qubit. Two 
current states of the three- Josephson-junction loop interact 
with the bias currents. Jo is the bias current applied across 
the capacitance C^, and fci's are the wave vectors of the 
Cooper pair wave function, (b) Plot for the effective potential 
t4ff({0±}) with / = 0.5 and (f)-=Q. The potential tilts due to 
a finite bias current /q. The sohd (dashed) line corresponds 
to positive (negative) Jq. 



with the constraint .g(0i) = 0. The Lagrange equation of 
motion with the above constraint produces the Kirchhoff 
equations in the qubit loop of Fig. [Ija): /o = /i + 12 and 

/l=/3. 

We introduce a coordinate transformation such as 
0+ = (02-f 03)/2 and 0- — (02 — 03)/2. In the usual fiux 
qubit experiment, the two Josephson junctions are nearly 
identical; thus, we set Ej2 = Ej^ — Ej. Although one 
can treat the general case numerically, this case gives a 
clear and intuitive picture for our qubit system. In this 
case, the effective potential is given by 



t^off({0±}) = -.BjiCos(27r/-20+)-2£;jcos( 
+^;^W -0+-0-). 



. COS( 



(11) 



Without the last term representing the linear coupling 
of the phase to the external bias current /q , the effective 
potential of Eq. (fTTj) can have local minima only for 
cos0_ = ±1, i.e., 0- = JTT, with j being an integer, 
where the qubit states are formed (See Fig. HJa)). 

Let the flux qubit be at the degeneracy point / = 0.5. 
When there is no bias current /q = 0, 0+ sa 7r/3(— vr/S) 
for the counterclockwise (clockwise) current state | J, 
) (I t)) with 0_ = 0, and the qubit energy levels cor- 
responding to the local minima of t/cff({0±}) are de- 
generate, E^ = E^. In this case, the qubit is opti- 
mally biased with respect to both the current /q and 
the magnetic flux /. A finite bias current tilts the effec- 
tive potential as shown in Fig. [2l^b), where the direction 
of energy level tilt depends on the sign of /q. Conse- 
quently, the effective potential of Eq. dTTI) is expressed 
as {E^ - $o/o«/27r)| i){i\ + {E^ + $o/oa/2^)| t)(t I 




FIG. 2: (Color online) (a) Effective potential of the current- 
biased flux qubit in the plane of (</>+, (f>-), where the effective 
potential of the local minima decreases as <j)+ or 0- increases. 
Here, we set 7o//c=0.05, / = 0.5 and Eji/Ej=0.8. (b) Plot 
for UeH{{4>±}) along the dashed line in the upper panel. 



In the basis of {|0), |1)}, the Hamiltonian is represented 
as 

huQ 



Ti. = —--(Tz + gsvnutax 



where hiijQ — 2t„ 



(16) 
and / 



uq is the qubit frequency, and a^^z 
are the Pauli and the identity matrices, respectively. The 
coupling strength g between the oscillating current and 
the qubit is given by 

5=^^a. (17) 

Near resonance, uj ss ojo, with a weak coupling g/h <C wq, 
one can apply the rotating wave approximation.^^ Then, 
we can observe a Rabi oscillation between the two qubit 
states |0) and|l) with a Rabi frequency i7p — g/h. 

Although we do not present it explicitly, in general the 
number of Josephson junctions can be any odd integer 
larger than one. In this case, at the two sides of the qubit 
loop divided by the bias current line, the numbers of 
Josephson junctions should be different from each other 
so that the symmetry of the flux qubit loop is broken. 
Then, the resulting coupling between the current and 
the phase of junction enables current-driven operation of 
the flux qubit. For the dc-SQUID (2-junction) qubit, the 
bias current flows across the Josephson junctions, but the 
symmetry of the dc-SQUID does not allow bias-current 
operation of the qubit states. 



apart from the constant, where 

Q« 10+1 w 7r/3. 



(12) 



The transitions between the two states, | J,) and | f); 
are induced by the charging energy given by the first 
term of the Lagrangian in Eq. ([5|). Using the tight- 
binding approximation, we can write the Hamiltonian for 
our qubit system as 






27r 



-«(u>ai-it>(ti), 



t)ai) 

(13) 



where tq is the transition rate between | 1) and | f) states. 
In order to operate the single qubit states, we apply a 
oscillating current. 



/o 



-/h sinwt. 



(14) 



As shown in Fig. [TJb), the effective potential vibrates 
with the frequency w, which produces oscillating diago- 
nal terms in Eq. (1131) . In transformed coordinates, these 
terms appear off-diagonal, giving rise to a Rabi oscilla- 
tion. The qubit states, |0) and |1), are the symmetric 
and the antisymmetric superpositions of | \.) and | '\), 
such that 



|o) = (ID + It))/^, 
|i) = (I ;> - 1 1»/^/2. 



(15) 



III. CIRCUIT QED 

The present current-biasing scheme for the flux qubit 
is implemented in the circuit QED architecture. In Fig. 
[31 the qubit is coupled to the resonator by an ac current 
flowing through a capacitance. The Lagrangian Cr of the 
transmission line resonator consists of the charge density 
mode q{x,t) and the current density mode I(x,t): 



Cr = 






(18) 



where I and c represent the inductance and the capaci- 
tance per unit length, and L the length of the resonator. 
Introducing the variable 6{x,t) = J_r,^dx'q{x',t), the 
Lagrangian becomes 
rL/2 
Cr= dx[il/2)idt9f - {l/2c){dx9f], (19) 

J-L/2 

where the voltage and the current of the resonator are 
given by V{x,t) = ^^^^ and I{x,t) = ^^^, respec- 
tively. For example, the voltage and the current of the 
resonator for the second harmonic mode can be repre- 
sented in terms of the boson creation and annihilation 
operators a^{t) and a{t) as^ 



-rrr \ tkJr f2TrX 



Ht) + a^{t)], 



I{x,t) = —i 



/ hcUr . ( ^TIX 

= —i\\ sm 

IL \ L 



[a{t)-aHt)], (20) 




FIG. 3: (Color online) Schematic diagram and equivalent 
lumped circuit representation of a three-j unction flux qubit 
coupled to a transmission line resonator. The flux qubit in- 
teracts with the current mode of the resonator. The arrows 
in the schematic diagram and equivalent circuit show the flow 
of oscillating current. As shown in the equivalent circuit, the 
oscillating current flows between the resonator and the qubit 
through the capacitance. Here, L is the length of the res- 
onator, and d is the width of qubit loop. 



where LOr = 2Trv/L and v — y^l/lc. The current profile 
for the second harmonic mode is shown in Fig. [3] 

In the circuit QED scheme with superconducting 
charge qubit, the qubit interacts with the voltage mode 
V{x,t). In contrast, the present current-biased flux 
qubit is coupled to the temporal fluctuation of local 
charge density, /(x,i), corresponding to the bias cur- 
rent applied to the qubit. In Fig. [3] the resonator and 
the qubit are coupled by a capacitance through the re- 
gion — (i/2 < X < d/2, where the charge fluctuation 
q{x, t) in the resonator produces the current flow into 

the qubit, /o(i) = J_^,^q{x,t)dx. As a result, follow- 
ing the current conservation condition in the resonator, 
q{x,t) = dl{x,t)/dx, and Eq. (I^Ul) . the current flowing 
into the qubit with width d is given by 



/o(i) = I{d/2,t)-I{-d/2,t) 



= -2i 



— sm -)[a(t) 



f(t)]. (21) 



The interaction between the flux qubit and the bias cur- 
rent Io{t) of the resonator is given by the last term of Eq. 
(HH). Then, from Eqs. dTS]) and ^, the total Hamiho- 
nian at the degeneracy point (E^ = E^) is given by 



H = fujJr ( CL a -\- 



1 



+ig{\\.){i 



-i,(u>(ti+it)ai) 

t)(t|)(a-at), (22) 



where the first term comes from the oscillating mode in 
the resonator. The first two terms of Eq. P^ disappear 



at the degeneracy point. In the basis of |0) and |1), the 
Hamiltonian is written as 



Ti, = hwrd a H — '^^z + ig<^x{o- ~ a )j 



(23) 



with hhJa = 2tq. Here, the single qubit gate is performed 
by a resonant external driving microwave fleld^. 

Since usually d ^ L, we have the expression of the 
coupling g as follows: 



g « q;<I>o 



huJr 

1l 



(24) 



Using the usual experimental values for the parameters^ 
such that d ~ i /im, L—b mm, lL—2.5 nH, ujr/2'K —15 
GHz, and a = tt/B, we estimate the coupling strength 
to be g/h ^ 120 MHz. For inductive coupling between 
the qubit and the resonator, the coupling strength, gic = 
^Ic sin a with $ being the magnetic flux threading the 
qubit due to the resonator current, can also be estimated 
to be 



5ic 



fj,Qd hw. 



2tiR V IL 



. sma, 



(25) 



where R is the mean distance between the resonator and 
the qubit loop and Ic ~ 2'kEj /<^o. With the same param- 
eter values, we see that 51c is smaller than g by an order 
of magnitude. We obtain the inductive coupling strength 
gic/h ~ 12 MHz with i? = 5 /zm and Ej/h = 200 GHz. 
Note that since the qubit is located at the nodal point 
of the current oscillation for the second harmonic mode, 
the inductive coupling at this point is negligible. The 
above value of inductive coupling has been obtained for 
the first harmonic mode. We have also assumed that 
the capacitance density is nearly uniform along the res- 
onator. In reality, the capacitance around the center of 
the resonator can be much higher due to the presence 
of the qubit loop, which can potentially provide a much 
stronger coupling g in the present scheme. 



IV. DECOHERENCE PROPERTY 

The qubit state of the present current-biased flux qubit 
is driven in a different way from that of the flux-driven 
flux qubit; thus, the decoherence property will also differ 
from each other. According to a recent review for the 
phase qubiti^, the dephasing times of the phase qubit and 
the three- Josephson junction flux qubit are comparable 
with each other. Since the only difference between the 
present current-biased flux qubit and the usual current- 
biased dc-SQUID qubit (phase qubit) is the number of 
Josephson junctions in the qubit loop, the decoherence 
property of both qubits can be analyzed in a similar way. 

In a recent experimenljii for the phase qubit, the ca- 
pacitance of Josephson junction is ~ 50fF while the three- 
Josephson junction flux qubit has a typically smaller 
value of capacitance of ^ 3fF ^^ with a small area of the 



Josephson junction. For both qubits, a large shunted ca- 
pacitance may reduce the decoherence from charge fluc- 
tuation. We employ the typical parameters of the three- 
Josephson junction flux qubit for the present current- 
biased qubit. Since both the capacitance and the critical 
current of the Josephson junction scale as the area of the 
junction, the critical current Ic of junction is also much 
smaller. 

For the present current-biased flux qubit, the dephas- 
ing due to the bias-current noise may cause decoherence 
of the qubit state, as in the phase qubiti^. The current 
noise is related to the 1// charge noise with a spectral 
density S'*(lHz)//. The phase noise is given by (0^) « 
[S'*(lHz)/CAC/](2/3) ln(1.778wioi)r^ where C is the ca- 
pacitance of the Josephson junction, A[/ is the barrier 
height of the potential, and wio is the qubit level spac- 
ing. Although for the present qubit, the spectral density 
S'*(lHz) and the capacitance C are small, these contri- 
butions cancel each other in (0^) because both the spec- 
tral density and the capacitance scale as the area of the 
junction^-. Further, the qubit level spacing wqi is similar 
between the two qubits. For the operation of the present 
qubit, we need not tilt the potential; thus, the barrier 
hight of the potential AC/ remains large. Hence, we con- 
sider the dephasing due to current noise in present qubit 
to still be comparable to the flux-driven flux qubit. On 
the other hand, the phase noise due to the critical-current 
fluctuation and to the flux 1// fluctuations is given by 
((/)2) « (5'|(lHz)L.7/At/)ln(0.401//™i)(woit)V6i2, with 
the Josephson inductance Lj = <i>o/27r/cCOS0 and the 
low- frequency cutoff fm- For the phase noise, the argu- 



ment is similar because both the spectral density S'KlHz) 
and the critical current Ic scale as the area of the junc- 
tion. 



V. SUMMARY 

We propose a new circuit QED architecture for the 
three- Josephson-junction flux qubit, where the flux qubit 
is coupled to the temporal charge density fluctuations of 
the transmission line resonator. The flux qubit is con- 
trolled by using a bias current. When three Josephson 
junctions are arranged asymmetrically, the energy levels 
of two qubit states with different chiralities couple to the 
external bias current. Rabi oscillations can be induced 
by an ac current at an optimal point with respect to both 
the bias current and the external magnetic flux. Remark- 
ably, the coupling between the qubit and the resonator 
is strongly enhanced compared to conventional inductive 
coupling. 
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